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Abstract. The impact of volcanic emissions, especially from 
passive degassing and minor explosions, is a source of uncer- 
tainty in estimations of aerosol indirect effects. Observations 
of the impact of volcanic aerosol on clouds contribute to our 
understanding of both present-day atmospheric properties 
and of the pre-industrial baseline necessary to assess aerosol 
radiative forcing. We present systematic measurements over 
several years at multiple active and inactive volcanic islands 
in regions of low present-day aerosol burden. The time- 
averaged indirect aerosol effects within 200 km downwind 
of island volcanoes are observed using Moderate Resolution 
Imaging Spectroradiometer (MODIS, 2002-2013) and Ad- 
vanced Along-Track Scanning Radiometer (AATSR, 2002- 
2008) data. Retrievals of aerosol and cloud properties at 
Kilauea (Hawai’i), Yasur (Vanuatu) and Piton de la Four- 
naise (la Reunion) are rotated about the volcanic vent to be 
parallel to wind direction, so that upwind and downwind re- 
trievals can be compared. The emissions from all three vol- 
canoes - including those from passive degassing, Strombo- 
lian activity and minor explosions - lead to measurably in- 
creased aerosol optical depth downwind of the active vent. 
Average cloud droplet effective radius is lower downwind 
of the volcano in all cases, with the peak difference ranging 
from 2-8 pm at the different volcanoes in different seasons. 
Estimations of the difference in Top of Atmosphere upward 
Short Wave flux upwind and downwind of the active volca- 
noes from NASA’s Clouds and the Earth’s Radiant Energy 
System (CERES) suggest a downwind elevation of between 
10 and 45 Wm“ 2 at distances of 150-400 km from the vol- 
cano, with much greater local (< 80 km) effects. Comparison 


of these observations with cloud properties at isolated islands 
without degassing or erupting volcanoes suggests that these 
patterns are not purely orographic in origin. Our observations 
of unpolluted, isolated marine settings may capture processes 
similar to those in the pre-industrial marine atmosphere. 


1 Background 

1.1 Aerosol indirect effects 

Aerosols affect the Earth’s albedo directly, through the ab- 

sorption and scattering of solar radiation, and indirectly, by 
altering the properties of clouds. Elevated levels of aerosol 
potentially lead to higher number densities of cloud con- 
densation nuclei (CCN), and for a parcel of cloudy air with 
a fixed mass of water, result in higher droplet concentra- 
tions and consequently smaller droplets and higher albedo 
(Twomey or first indirect effect; Twomey, 1977). In addi- 
tion, smaller cloud droplets may result in the suppression of 
precipitation and therefore longer cloud lifetime, i.e. higher 
albedo (second indirect effect, Albrecht, 1989), although the 
importance of this effect is thought to vary between atmo- 
spheric regimes (Stevens and Feingold, 2010). If the overly- 
ing air is sufficiently dry, the evaporation of small droplets in 
a polluted cloud is enhanced, so that cloud water content de- 
creases as droplet concentrations increase (Ackerman et al., 
2004). The presence of aerosol in a cloud may also cause 
the evaporation of cloud droplets due to aerosol absorption 
of solar radiation, in addition to several less well-understood 
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perturbations to droplet character in mixed phase clouds (e.g. 
Hansen et al., 1997). The impact of secondary aerosol indi- 
rect effects has not been well quantified or verified by obser- 
vations (Rosenfeld et al., 2014). 

Aerosol indirect effects (AIE) encompass the combined 
effects of both anthropogenic and natural emissions. The 
contributions of AIE to the Earth’s radiative balance are 
strong, yet also highly uncertain. Carslaw et al. (2013) sug- 
gest that 45 % of variance in post- 1750 aerosol forcing (the 
contribution of aerosols to the change in the Earth’s radiative 
balance as CO 2 levels rise) is from natural sources, which 
are hard to isolate and measure in the polluted present-day 
atmosphere (Andreae , 2007). The response of cloud prop- 
erties to additional CCN is non-linear: excess CCN has the 
greatest effect on albedo when background levels of aerosol 
are lowest (e.g. Lohmann et al., 2005). Present-day AIE are 
therefore expected to most closely match pre-industrial pro- 
cesses in remote marine environments where pollution levels 
are low (Andreae , 2007). 

Direct observations of the local effects of aerosols on 
clouds have so far been dominated by measurements of ship 
tracks (e.g. Durkee et al., 2000; Ackerman et al., 2000; 
Campmany et al., 2009; Christensen and Stephens, 2011). 
Ship tracks provide an ideal experiment for isolating the im- 
pact of aerosols, as the polluted clouds are similar in origin 
and thickness to clean clouds in the surrounding cloud deck, 
although they may be elevated in height (Christensen and 
Stephens, 2011). The occurrence and albedo of ship tracks 
has been shown to depend on pre-existing cloud structure, 
height and humidity (e.g. Chen et al., 2012). However, Pe- 
ters et al. ( 2011 ) found no statistically significant impact of 
aerosol from shipping on a large scale, away from the ship 
tracks themselves. 

1.2 Volcanic aerosol indirect effect 

The impact of volcanogenic aerosol on cloud properties de- 
pends on injection height into the atmosphere, pre-existing 
aerosol burden and synoptic conditions. Explosive eruptions 
can inject aerosol directly into the stratosphere (as coarse par- 
ticles from fragmented magma, vent wall erosion and con- 
densation of magmatic gases) and provide the gas-phase re- 
actants for the production of finer aerosols in the plume or 
ambient atmosphere (e.g. sulfate aerosol). Direct and indirect 
aerosol effects associated with explosive eruption products 
in the stratosphere have been relatively well characterised 
(Sassen, 1992; Robock, 2002; Schmidt et al., 2010). How- 
ever, the role of aerosol in the troposphere from “passive” de- 
gassing and minor eruptions is less well understood (Mather 
et al., 2003; Oppenheimer et al., 201 1). Time-averaged emis- 
sions from passive degassing are thought to make up a 
high proportion (~ 30-70%) of the volcanic SO 2 flux to 
the atmosphere (Andres and Kasgnoc, 1998; Halmer et al., 
2002; Mather et al., 2003). Volcanogenic aerosols in the tro- 
posphere are expected to be dominated by sulfate aerosol 


(SC> 4 _ ) produced from the oxidation of SO 2 by OH - , H 2 O 2 
or O 3 above the boundary layer, or in high-temperature re- 
actions of SO 2 or H 2 S in the volcanic vent (e.g. Mather 
et al., 2006). Ambient-temperature SO 2 reactions can take 
place over several days, so the impact of volcanic aerosol on 
number density of CCN may extend hundreds of kilometres 
away from the aerosol source (Eatough et al., 1995). Volcanic 
emissions are one source of natural aerosol that, in combina- 
tion with anthropogenic aerosol, result in AIE. For readabil- 
ity, we refer to this contribution as “volcanic aerosol indirect 
effects” (VAIE). 

The flux of volcanic aerosol to the present-day troposphere 
is lower than the anthropogenic flux, but its impact on the 
radiative budget may be disproportionately significant, as 
volcanic gases are commonly emitted at heights above the 
boundary layer into the free troposphere so that aerosol life- 
times may be longer (e.g. Graf et al., 1998). In global model 
simulations of present-day CCN, uncertainties in volcanic 
aerosol are lower than those for biomass burning and anthro- 
pogenic SC> 4 “ sources (e.g. Lee et al., 2013). However, vol- 
canic emissions contribute a lower proportion of CCN to the 
atmosphere today than in pre-industrial times, so the quan- 
tification of VAIE are particularly important for estimating 
the baseline cloud radiative state (Schmidt et al., 2012). Vol- 
canic SO 2 emissions have the greatest potential range of all 
the natural emissions required to estimate global mean forc- 
ing uncertainty (Carslaw et al., 2013). 

Ground-based measurements of plumes near passively 
degassing vents have captured (1) an abundance of SO 4 - 
aerosol larger than the critical diameter to act as a CCN at 
typical supersaturations (e.g. Allen et al., 2002; Martin et al., 
2008; Mather et al., 2012) and (2) particle growth attributed 
to the condensation of water vapour onto sulfate particles be- 
tween 0 and 20 km distance from the vent (Villarica, Chile; 
Mather et al., 2004). Satellite remote sensing studies have 
gone a step further to measuring VAIE during periods of 
high flux degassing. Gasso (2008) shows alteration to stra- 
tocumulus properties (“volcano tracks”) in Moderate Reso- 
lution Imaging Spectroradiometer (MODIS) and Advanced 
Microwave Scanning Radiometer - Earth Observing System 
(AMSR-E) images where large plumes from volcanoes in 
the South Sandwich islands and the Aleutians interact with 
marine boundary layer clouds. Instances of increased cloud 
brightness in the presence of an obvious volcanic plume were 
identified by browsing MODIS visible images from 6 months 
in 2006, and cloud property data from a few selected days 
were analysed to demonstrate first and second AIE in the 
volcano tracks. Yuan et al. (201 1) have also assessed the im- 
pact of high degassing flux associated with a new vent open- 
ing at Halema’uma’u crater, Hawai’i, in 2008 on trade cu- 
mulus clouds. The authors compare retrievals of cloud and 
aerosol properties averaged over 3 months from inside the 
Halema’uma’u plume to properties both outside the plume 
over the same time, and to the average properties retrieved 
over a much longer period. Both studies (Gasso, 2008; Yuan 
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et al., 2011) find evidence for increased albedo and for in- 
creased cloud lifetime in marine clouds. 

1.3 Motivation and aim 

Although volcanic aerosol acting as CCN has been observed 
at several volcanoes (e.g., Mather et al., 2004; Gasso, 2008; 
Yuan et al., 2011), its net effect has not previously been quan- 
tified over extended periods of time. More representative 
measurements are important for testing and verifying the pre- 
dictions of global microphysical aerosol models (e.g. Martin 
et al., 2006; Dentener et al., 2006; Schmidt et al., 2012) that 
rely on inventories of sulfur emission such as the one com- 
piled by Andres and Kasgnoc (1998). 

As measurements of VAIE so far have been made over 
short time intervals (3-5 months) and during elevated activ- 
ity - e.g. frequent MODIS thermal anomalies at Montagu 
and Saunders Islands (Gasso, 2008) and Halema’uma’u vent 
opening (Yuan et al., 2011) - they are not likely to be rep- 
resentative of the long-term impact of “background” activity 
on tropospheric cloud properties. 

This study presents an approach for detecting VAIE for the 
particular case of isolated, active volcanic islands. We use 
averages of 6-10 years of data and make no selection on the 
basis of activity, so that our results may be considered rep- 
resentative of the net effect of volcanic emissions on clouds 
at a particular volcano. Three “control” islands (no volcanic 
aerosol) are also considered, to provide a comparison with 
possible island (e.g. orographic) effects on cloud properties. 
Our aim is to capture the VAIE associated with “background 
activity” - that is, periods when degassing or other persistent 
activity results in aerosol emission into the troposphere and 
that may not result in visibly identifiable “volcano tracks”. 


2 Satellite data 

Our analysis uses satellite-retrieved aerosol and cloud prop- 
erties, specifically aerosol optical depth (AOD) at 550 nm, 
cloud mid- visible optical depth (COD), and cloud droplet ef- 
fective radius (CER). Two satellite products are used for each 
of these quantities. The first is the Collection 6 MODIS At- 
mospheres Level 2 Joint product (MYDATML2), from the 
Aqua satellite (data from 2002-2013). These aerosol/cloud 
data products are described by Platnick et al. (2003), Acker- 
man et al. (2008) and Levy et al. (2013). MODIS cloud re- 
trievals (resolution 5x5 km 2 ) are included where the cloud 
fraction from MODIS cloud product is >0.2, and Qual- 
ity Assurance (QA) values are >0 (removing data where 
confidence in the retrieval was low). For aerosol retrievals 
(resolution 10 x 10 km 2 ), we also require QA >0. The sec- 
ond aerosol/cloud data source is the Oxford-RAL Aerosols 
and Clouds (ORAC) Global Retrieval of Cloud Parameters 
and Evaluation (GRAPE) data set, from Advanced Along- 
Track Scanning Radiometer (AATSR) measurements, from 


2002-2008 (Thomas et al., 2009; Sayer et al., 2010, 2011; 
Poulsen et al., 2012), with data quality-filtered as recom- 
mended by the above. Some significant differences in CER 
between these two data sets were observed by Sayer et al. 
(2011), attributed in part to the different wavelengths used 
by the retrieval algorithms, and consequent differences in 
sensitivity to cloud vertical structure (e.g. Platnick, 2000). 
For both aerosol data sets, we require pixel cloud fraction 
<0.8. Although there may be artefacts in AOD retrieved in 
conditions of broken cloud (e.g. Quaas et al., 2010; Grandey 
et al., 2013), we use AOD retrievals when cloud fraction is 
up to 0.8 to maximise the number of retrievals in our analy- 
sis. Cloud fraction, cloud top pressure and temperature were 
used to distinguish between retrievals from different atmo- 
spheric heights and under different synoptic conditions. We 
restrict our observations to liquid water rather than ice clouds 
(and to data where cloud top pressures were >440 mb). All 
properties are resampled to a grid of 10 km resolution to 
simplify the comparison. Additionally, data were stratified 
by season to account for potential large-scale differences in 
aerosol/cloud properties and data sampling rates in different 
seasons. 

We use Clouds and the Earth’s Radiant Energy System 
(CERES) Single Scanner Footprint data (Top of Atmosphere 
upward Short Wave flux, resolution 20 x 20 km 2 ; Geier et al., 
2003) to estimate radiative impact separately from these mi- 
crophysical property retrievals. Although differences in solar 
zenith (expected to be < 10° in all cases and < 5° for the ma- 
jority of retrievals) may result in variations in SW flux across 
the region surrounding the volcano, correlations with either 
cloud properties or wind direction are improbable and are 
therefore unlikely to affect comparison of upwind and down- 
wind SW fluxes. 

Data are selected within a square of side length 4° (up 
to 600 km, depending on latitude) centred on the volcanic 
vent. We then convert the latitude and longitude for each re- 
trieval pixel into polar coordinates with the origin located 
on the volcano’s summit, as recorded in the Smithsonian 
database (Siebert and Simkin, 2002). Horizontal wind ve- 
locity components from the European Centre for Medium- 
Range Weather Forecasts (ECMWF, ERA-Interim, spatial 
resolution =1.5°, e.g. Dee et al., 2011) are used to rotate 
the aerosol and cloud properties according to wind direc- 
tion, so that they can be plotted according to their expected 
position upwind or downwind relative to the volcano. ERA- 
Interim horizontal wind velocity components (available five 
times a day) are selected for the time of day closest to the 
satellite overflight. We assume that aerosols and their precur- 
sor gases are emitted from each volcano at approximately its 
summit height, so wind data are selected for a pressure level 
expected to lie just above this height. The volcano’s summit 
or the most active vent is treated as the origin for rotation 
because for all cases here, gases were emitted from multiple 
vents contributing varying proportions of total aerosol emis- 
sions. 


www.atmos-chem-phys.net/14/10601/2014/ 


Atmos. Chem. Phys., 14, 10601-10618, 2014 


10604 


S. K. Ebmeier et al.: Passive volcanic degassing and cloud properties 


The resulting averaged images of rotated cloud and aerosol 
properties allow us to examine systematic trends that are not 
apparent in the properties retrieved for individual days, or in 
the average of all our data before rotation. Unlike previous 
studies at periods of high SO 2 emission (e.g. Gasso, 2008; 
Yuan et al., 2011), volcano tracks were rarely identifiable in 
individual MODIS visible images or in cloud properties on 
specific dates. For MODIS and AATSR, retrievals of cloud 
and aerosol properties are mutually exclusive, as aerosols are 
retrieved only in clear-sky conditions. CERES Top of Atmo- 
sphere (ToA) upward Short Wave (SW) flux is available for 
both cloudy and clear-sky conditions. 

2.1 Choice of targets 

We select target volcanoes where retrieval uncertainties are 
expected to be low, to allow us to identify low-magnitude 
perturbations in cloud properties. Uncertainties in AOD re- 
trieved from MODIS data over land are thought to be on av- 
erage three times greater than over water, where models of 
surface reflectance are better (Remer et al., 2005). We there- 
fore focus our study on isolated volcanic islands, avoiding 
the higher retrieval uncertainties, greater variability in cloud 
characteristics associated with continents and the systematic 
dependence of cloud form on wind over coastlines or high 
topography (e.g. Brenguier et al., 2003). The ideal target is 
a small and low-lying island volcano, with a high, persistent 
degassing flux and located a long way from other sources of 
tropospheric aerosol, such as cities, heavy industry or ship- 
ping. 

Three isolated, persistently if variably active volca- 
noes were chosen to test our approach (Figs. 1 and 2): 
Kilauea (Hawai’i), Yasur (Vanuatu) and Piton de la Four- 
naise (Reunion). As described above, measurements of a 
VAIE downwind of Kilauea during the opening of the new 
Halema’uma’u vent were previously made by Yuan et al. 
(2011). Although the results presented here are also for 
mafic, persistently active volcanoes, they capture a range of 
eruption styles and degassing fluxes (Table 1). 

A VAIE is expected to be characterised by (a) elevated 
AOD and (b) depressed CER downwind of the volcanic 
source relative to upwind. However, we also expect an oro- 
graphic effect over isolated peaks in topography (e.g. Jiusto, 
1967) to cause a systematic upwind/downwind difference 
in CER: when moist air is lifted up over an island it ex- 
pands, cools and condenses, potentially causing precipitation 
on the windward side of the island. We therefore also con- 
sider data for three “control” islands known not to emit vol- 
canogenic aerosol or precursor gases (Tristan da Cunha in 
the South Atlantic; Ofu-olosega, American Samoa and Fiji, 
both in the South Pacific - Table 1), chosen for their similar- 
ities in height and diameter to our target volcanoes (Fig. 2), 
and their isolation from other sources of volcanic or anthro- 
pogenic aerosol. 


3 Results 

The differences between upwind and downwind aerosol and 
cloud properties are much greater at the volcanoes than at 
the control sites. AOD is elevated, CER suppressed and ToA 
upward SW flux higher downwind of all three volcanoes 
(Figs. 3-5; Table 2). None of the “control” islands show a 
notable difference between upwind and downwind properties 
or a difference in SW flux exceeding ±10 Wm" 2 (Fig. 5). 

Degassing processes and patterns are notably different at 
each of the volcanoes, and we discuss results from each vol- 
cano below with reference to other measurements of volcanic 
aerosol and prevailing meteorological conditions (Sects. 3.1- 
3.3). 

3.1 Kilauea, Hawai’i 

Kilauea is the most active of the Hawaiian volcanoes and 
contributes 4 % of the total emissions in Andres and Kasg- 
noc’s (1998) time-averaged compilation of emissions from 
49 continuously erupting volcanoes. Activity is effusive, and 
split between the East Rift zone, where lava effusion has 
been semi-continuous since 1983, and periodic degassing, 
lava pond activity and occasional small explosions at the vol- 
cano’s summit. Prior to 2008 approximately 90% of SO 2 
emissions were from the East Rift zone, and between 2002 
and the beginning of 2008 average total annual emissions 
were 6.58 x 10 5 Mgyr -1 (Elias and Sutton, 2012). A com- 
parison of the USGS ground-based spectrometer measure- 
ments with MODIS AOD prior to 2008 for days where both 
measurements were possible (61 days) does not show a re- 
lationship between them (Fig. 6), perhaps due to the time 
difference between the measurements, the lag between SO 2 
emission and aerosol formation or processes that affect the 
concentrations of aerosols, such as dilution. 

The opening of a new vent at Halemau’ma’u crater in 
2008 was accompanied by a rise in total SO 2 flux to 
7000 Mg day -1 , constituting a doubling of annual SO 2 emis- 
sion rate from 2007 to 2008. Since 2008 summit emissions 
steadily increased as a proportion of total SO 2 emission 
from Kilauea, though total SO 2 emission had dropped to 
~ 2000 Mg day -1 by 2010 (Elias and Sutton, 2012). Yuan 
et al. (2011)’s measurements of the Halema’uma’u 2008 
plume demonstrate a significant increase in cloud fraction 
and perturbation to cloud properties (a) inside the aerosol 
plume, relative to outside and (b) during elevated emis- 
sion in 2008, relative to the mean values for the 11 -year 
Terra archive. On 17 July 2008 Yuan et al. (2011) found 
an average background value of CER of 1 8 pm, reduced 
to ~ 13 pm within the aerosol plume 800-1000 km down- 
wind of Kilauea. We find a similar magnitude of reduc- 
tion in CER over a much smaller area (just ^70 km down- 
wind) in averaged MODIS data from 2002-2007, when de- 
gassing fluxes were lower, as well as in the multiannual aver- 
age 2002-2013 (Table 2, Fig. 3). Yuan et al. (2011) found 
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MODIS 550 nm AOD composite 



Figure 1. Multiannual mean (2002-2013) of Collection 6 MODIS Aqua AOD at 550 nm (merged data set). Grey indicates that there are no 
data. Volcano locations are indicated by blue diamonds, “control” islands by blue crosses. 
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Figure 2. (a-f) Dimensions, topography and seasonal distribution of data sampling for each volcano and control island in the study. Wind 
roses show daily ECMWF wind directions, black bars indicate the directions towards which the wind was blowing. Locations of degassing 
vents are marked by red stars. Bar charts show the number of retrievals used in this study for each month of the year, Squares represent 
number of ENVISAT retrievals, circles, MODIS. Filled shapes are cloud and outlines are aerosol. 


a 20 Wm -2 perturbation within the 2008 Halema’uma’u 
plume. At distances of ~ 150-400 km downwind of Kilauea 
we find perturbations of between 20 and 45 Wm ~ 2 (mean 
= 28Wm~ 2 ) between 2002 and 2013, well constrained to 
a region within ~45° of the downwind direction. Within 
^80 km of the vent, where many retrievals would have been 
over Hawai’i island itself, the perturbation is even higher, ex- 
ceeding 90 Wm -2 (Fig. 5a-b). 


Sulfates dominate the compositions of aerosols measured 
at Kilauea’s summit (modal diameter 0.44 pm), and make up 
~ 1 % by mass of plume SO 2 concentrations (Mather et al., 
2012). Porter et al. (2002) made sun photometer and Li- 
DAR measurements of SO 4 - 9 km downwind of a vent on 
the East Rift zone in 2001 and found an aerosol dry mass 
flux rate of 53 Mg day -1 . They inferred a half-life of 6 hours 
for SO 2 in the plume. Assuming an average wind speed of 
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Table 1 . Location, physical properties, degassing flux and eruptive character for volcanoes and “control” islands discussed in this study. 


Volcano/Island 

Fatitude 

(°) 

Longitude 

(°) 

Maximum diameter 
of island (km) 

Emission (summit) 
height (m) 

“Background” SO 2 
flux* (Mg day -1 ) 

Activity 2002-2008 

Volcanoes 

Kilauea 

19.42 

-155.29 

130 

1222 (4170) 

1800 (2002-2007) 1 

Effusive lava flows, lava lake 

Yasur 

-19.53 

169.44 

40 

361 

633 (2004-2008) 2 

Strombolian 

Piton de la Foumaise 

-21.23 

55.71 

70 

2632 

Inter-eruptive flux is very low 3 

Minor explosions and basaltic flows 4 

“Control Islands 44 

Fiji 

-17.63 

178.02 

140 

1324 

- 

- 

Ofu-olosega 

-14.18 

169.62 

9 

639 

- 

- 

Tristan da Cunha 

-37.09 

-12.28 

12 

2060 

- 

- 


* “SO 2 flux” refers to time-averaged values of emission not associated with explosions as found from measurements presented by 1 Elias and Sutton (2012), 2 Bani et al. (2012) 
and 3 Coppola et al. (2009); Garofalo et al. (2009) and Di Muro et al. (2012). 

4 Periods when eruptions at Piton de la Foumaise occurred: 10.2010-12.2010, 1 1.2009-01.2010, 09.2008-02.2009, 07.2006-05.2007, 02.2005, 05.2004-10.2004, 
05.2003-01.2004 and 11.2002-12.2002; Siebert et al. (2010). 


range 5-10 m s -1 , if correct and broadly applicable in differ- 
ent atmospheric conditions, this implies that SO 2 mass falls 
below 10 % of its original value by a distance of 90-180 km 
downwind. Our estimations of average AOD place peak val- 
ues at ~ 100 km downwind of Kilauea’s summit, consistent 
with oxidation of SO 2 as the primary source of aerosol. 

Aerosol plume dispersal in Hawai’i is dominated by rel- 
atively stable trade winds from the northeast. Clouds are 
mostly trade cumuli, which are typically capped by the trade 
wind inversion at a height of a few kilometres. Although 
Kilauea’s summit elevation is only 1222 m, it overlaps the 
lower flanks of the much larger, but less active, Mauna 
Loa (4170 m). There is a notable difference in the num- 
ber of MODIS retrievals flagged as cloud downwind rela- 
tive to upwind (retrievals with cloud fraction > 0.2 are shown 
in Fig. 7). More cloud retrievals were made downwind of 
Kilauea, with the difference peaking at a distance of about 
60 km downwind of the vent. Elevated downwind COD is 
likely to be associated with both the formation of orographic 
cloud and the action of volcanogenic aerosol as CCN. 

3.2 Yasur, Vanuatu 

Yasur is the southernmost of the Vanuatu island arc’s ac- 
tive volcanoes and exhibits almost continuous Strombolian 
to Vulcanian activity in an eruption that has lasted at least 
300 years. Magmatic gases, mostly SO 2 , are continuously re- 
leased from three vents in Yasur’s crater and typically rise to 
heights of 700-900 m before being carried to the northwest 
by trade winds (Bani et al., 2012). Measurements made by 
Bani and Lardy (2007) and Bani et al. (2012) on 33 different 
days between 2004 and 2008 find an average degassing flux 
of 633 Mg day -1 SO 2 (standard deviation of daily measure- 
ments = 270 Mg day -1 ). Bani and Lardy (2007) estimate that 
in 2004-2005 Yasur contributed 1-2 % of estimated global 
time-averaged volcanic SO 2 emissions to the troposphere, 
close to Andres and Kasgnoc (1998)’s estimation of 3 % pre- 
1990s. 


The nearest other aerosol sources to Yasur are just un- 
der 400 km away, at the volcanoes of Lopevi (156 Mg day - 1 
SO 2 ) and Ambrym (5440 Mg day -1 SO 2 ) (Bani et al., 2012). 
Yasur is close to being a sea-level point SO 2 source, with a 
summit elevation of just 361 m. We therefore expect any oro- 
graphic effects to be minimal. There are also no indications 
of orographic clouds in rotated cloud properties at the small, 
remote islands in American Samoa (Ofu-olosega, elevation 
639 m) used as “control” sites. 

We measure elevated AOD (e.g. seasonal average upwind- 
downwind difference of 0.03-0.06 for MODIS (Aqua); 0.02- 
0.03 for AATSR) downwind of Yasur, relative to that upwind. 
Ground-based SO 2 measurements at Yasur have been too in- 
frequent (Bani et al., 2012) to allow a useful comparison with 
AOD retrieved from MODIS data. There is also a net differ- 
ence in CER between the downwind and upwind sectors (sea- 
sonal average differences are —3 — 4 pm, MODIS; 0 — 4 jam 
AATSR). Although some of this difference can be attributed 
to the impact of excess aerosol from Yasur, this is superim- 
posed on a regional trend in aerosol and does not clearly in- 
dicate statistically significant AIE. Similarly, the elevation 
in ToA upward SW flux downwind of Yasur (Fig. 5c-d) is 
likely to be influenced by regional variation in cloud proper- 
ties, rather than purely volcanic effects. 

3.3 Piton de la Fournaise, Reunion 

Piton de la Foumaise, Reunion, differs from the previous 
two examples in that intra-emptive SO 2 flux is low (Khokhar 
et al., 2005; Di Muro et al., 2012). However, Piton de la Four- 
naise is very active, on average it empts or experiences an 
intmsion once every 8 months (Peltier et al., 2009). Erup- 
tions are sometimes associated with measurable SO 2 emis- 
sion (Khokhar et al., 2005; Bhugwant et al., 2009), and emis- 
sion is of longer duration during periods of distal lava effu- 
sion when lava entering the sea also results in a high flux 
of water vapour to the atmosphere (e.g. Gouhier and Cop- 
pola, 2011). Between emptions, degassing is primarily from 
low temperature fumeroles and fractures around the active 
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Kilauea Yasur Piton de la Foumaise Tristan da Cunha America Samoa Fiji 




Figure 3. Profiles of average aerosol and cloud properties upwind and downwind of the degassing volcano: Kilauea, Hawai’i; Yasur, Vanuatu 
and Piton de la Foumaise, Reunion and at “control” islands: Tristan da Cunha, South Atlantic; American Samoa, West Pacific and Fiji, West 
Pacific. Profile values are the average of arcs of j at equal distance upwind and downwind from the volcano vent for MODIS Aqua data from 
between 2002 and 2013. Separate profiles for different time of year are shown in different colours: December-Febmary in red, March-May 
in blue, June-August in green and September-November in black. Error bars show one standard error for this average and are generally larger 
for low sample sizes. Positive values on the v-axis are downwind, negative are upwind. The maximum diameter of the islands is indicated 
on the parameter plots as black dashed lines. 
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Figure 4. Profiles of average aerosol and cloud properties from AATSR between 2002 and 2008, upwind and downwind of the degassing 
volcano, presented as in Fig. 3, above. 
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summit crater (Di Muro et al., 2012). Again, unlike Kilauea 
and Yasur, these weak emissions are dominated by water 
vapour, and hydrogen sulfide (H 2 S) is the main sulfur species 
(Di Muro et al., 2012). La Reunion is located near an impor- 
tant shipping lane (Peters et al., 2011), so it is possible that 
some anthropogenic sources of sulfate aerosol may also be 
present. 

Piton de la Foumaise is the only one of the volcanoes in- 
vestigated here that is likely to have had a high SO 2 flux asso- 
ciated with explosive eruption. The SO 2 released during the 
collapse of the Dolomieu crater in April 2007 (estimated to 
be 935 ±244 kilotons by Gouhier and Coppola, 2011) rose 
quickly to heights > 3000 m (above the trade winds), and 
reached distances of 800-1000 km away from La Reunion 
within days (Tulet and Villeneuve, 2011). Where SO 2 emis- 
sion accompanies explosions, and aerosol is carried above 
the trade winds, the aerosol and associated impact will be 
spread over a much greater area. At distances of ~ 1000 km, 
the plume from this event reached heights of around 10 000 m 
above sea level (Tulet and Villeneuve, 2011), so aerosol 
from the plume are only likely to settle out and re-enter the 
lower atmosphere at great distance from the volcano. Piton 
de la Fournaise’s other six eruptions between 2002 and 2008 
(Siebert et al., 2010) were all less explosive, emitting aerosol 
lower into the atmosphere. 

Differences between upwind and peak downwind AOD 
at Piton de la Foumaise range from 0.01 to 0.08 (MODIS) 
in different seasons at ~ 50 km distance from the volcano’s 
summit, while differences in CER are —2 — 10 pm (seasonal 
and multiannual mean MODIS and AATSR values shown in 
Figs. 3, 4 and 8.) Peak downwind elevation in To A upward 
SW flux is ~ 28 Wm -2 at 20 km downwind of the volcano 
(Fig. 5e-f). 

Our observations capture the impact primarily of low- 
altitude aerosol emission. However, our observations of 
VAIE are not limited to inter-emptive emission. In Fig. 9 
aerosol and cloud properties are split according to whether 
or not Piton de la Foumaise was empting (as recorded in the 
Smithsonian database; Siebert and Simkin, 2002). Aerosol is 
elevated all around the volcano during periods of emption 
(compare Fig. 9a and b), but particularly in the downwind 
quadrant. Similarly, CER is lowest during periods of emp- 
tion in the downwind sector (compare Figs. 9c and d). 

3.4 Distinguishing between orographic and volcanic 
effects 

Orographic clouds form when moist air is forced upwards by 
high topography so that it expands and cools below its satu- 
ration point. Their form depends on relative humidity, wind 
speed and the height and geometry of topography. At high 
humidity and lower wind speeds, droplets condense above 
the windward slopes of a mountain, sometimes causing oro- 
graphic precipitation. Wave clouds may also form in per- 
pendicular bands downwind of high topography. At isolated, 


steep-sided mountains high wind speeds can cause local up- 
lifting on the upper slopes in the lee side of the mountain 
(“banner cloud”). As the formation of orographic cloud is 
controlled by wind strength and direction, its presence is ex- 
pected to introduce systematic features to our rotated cloud 
properties. 

For a tall, isolated island, orographic clouds may form (1) 
over high land upwind of the summit and (2) at the crests 
of lee waves (rotor and lenticular clouds), downwind and ex- 
tending hundreds of kilometres away from the island (e.g. 
Houze, 1994). We observe elevated COD in averaged cloud 
property data over the islands with high topography (Kilauea, 
Reunion, Fiji and to some extent Tristan da Cunha). COD 
is elevated immediately upwind of the islands’ summit and 
remains high over land downwind. Any difference in COD 
between the windward and lee sides of the islands at greater 
distance is very low, suggesting that the contribution of wave 
clouds to the average values is small. 

This is confirmed by estimating the average cloud proper- 
ties for days with different tropospheric stabilities (Fig. 10). 
Tropospheric stability was approximated by the unsaturated 
moist Froude number estimated from ECMWF reanalysis 
data: Fr 2 = where u is wind speed at the height ( h ) 
of the island and N is the unsaturated moist Brunt-Vaisala 
frequency. N is also estimated from ECMWF temperatures 
at atmospheric pressures >750 mb (A 2 = |-^r, where g = 
gravitational acceleration, 0 V = virtual potential temperature 
and z = atmospheric thickness). Where Froude number is low 
(< 0.1, blue on Fig. 10) air will flow around a mountain, 
rather than over it. For Froude numbers closer to 1 (0. 1-0.5, 
green and >0.5 red on Fig. 10), lee waves and wave clouds 
will form. At higher Froude numbers air flows over the is- 
land very fast and there are no stable layers. Although COD 
is higher at greater Froude number at Yasur and Piton de 
la Foumaise (Fig. 10), atmospheric stability seems to have 
a limited effect on systematic differences between CER up- 
wind and downwind. 

Orographic cloud has smaller droplets at a higher den- 
sity than marine clouds. For example, Jiusto’s (1967) air- 
craft sampling of both marine and orographic cloud upwind 
of Kilauea found that marine clouds had an average radius 
of 2 2±3pm (droplet concentration 45 ±22 cm -3 ) while 
clouds on the flanks of Mauna Loa had an average radius of 
14 ± 3 pm (droplet concentration 100 ± 50 cm -3 ). Although 
the values are not directly comparable, we find a similar de- 
crease in droplet size (seasonal averages —2 — 8 pm, see Ta- 
ble 2) on the windward side of Big Island Hawai’i. However, 
CER does not reach its minimum value until almost 100 km 
downwind, and remains lower than upwind values until al- 
most 400 km from land (see Fig. 3). We see no such effect 
in clouds downwind of Fiji (also ~ 100 km across and in the 
path of easterly trade winds, but lower in elevation) or Amer- 
ican Samoa, but Tristan da Cunha shows a small decrease in 
droplet size in data from December-February. We note that 
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Table 2. Average difference in upwind and downwind cloud and aerosol properties for MODIS (Aqua) and AATSR retrievals for volcanoes 
discussed in this study. Uncertainties are the square root of the sum of squares of the standard errors of upwind and downwind values. 


Volcano 



A AOD 1 

A COD 2 

A CER (nm) 1 

Peak distance (km) 3 

Piton de la Foumaise 

MODIS 

DJF 

0.08 ±0.03 

2 ± 1 

— 8 ± 1 

-25 



MAM 

0.03 ±0.03 

3 ± 1 

— 3 ± 1 

-30 



JJA 

0.03 ±0.02 

2 ± 1 

— 3 ± 1 

-50 



SON 

0.04 ±0.02 

<1 

— 4 ± 1 

-40 


AATSR 

DJF 

0.02 ±0.02 

<1 

—4 ±2 

-40 



MAM 

0.03 ±0.02 

—2 ±4 

—3 ±2 

-15 



JJA 

0.03 ±0.02 

<1 

—2 ±2 

-25 



SON 

0.01 ±0.02 

<1 

—4 ±2 

- 

Yasur 

MODIS 

DJF 

0.06 ±0.01 

<1 

—4 ± 1 

-40 



MAM 

0.03 ±0.01 

2 ± 2 

— 3 ± 1 

-25 



JJA 

0.04 ±0.01 

5 ± 2 

— 4 ± 1 

-35 



SON 

0.05 ±0.01 

3 ±2 

— 3 ± 1 

-40 


AATSR 

DJF 

0.03 ±0.02 

<1 

<1 

- 



MAM 

0.02 ±0.02 

2 ± 3 

—4 ±2 

-40 



JJA 

0.02 ±0.02 

5 ± 3 

<1 

- 



SON 

0.03 ±0.02 

2 ± 3 

—2 ±2 

-40 

Kilauea 

MODIS 

DJF 

0.06 ±0.01 

2 ± 1 

— 6 ± 1 

-90 



MAM 

0.11 ±0.02 

4 ± 1 

— 6 ± 1 

-100 



JJA 

0.14 ±0.02 

4 ± 1 

— 8 ± 1 

-70 



SON 

0.05 ±0.02 

2 ± 1 

— 8 ± 1 

-90 


AATSR 

DJF 

0.07 ±0.03 

4 ± 2 

— 2 ± 1 

-100 



MAM 

0.08 ±0.03 

5 ± 3 

-4 ±3 

-90 



JJA 

0.08 ±0.02 

5 ± 3 

— 5 ± 3 

-90 



SON 

0.06 ±0.03 

2 ± 2 

—2 ±2 

-90 

Tristan 

MODIS 

DJF 

-0.02 ±0.02 

2 ± 1 

— 2 ± 1 

-20 



MAM 

0.01 ±0.02 

3 ± 1 

-1±1 

-20 



JJA 

-0.01 ±0.02 

3 ± 1 

— 2 ± 1 

-20 



SON 

0.03 ±0.02 

2 ± 1 

-1±1 

-20 


AATSR 

DJF 

-0.01 ±0.03 

1 ±2 

—2 ±2 

- 



MAM 

-0.01 ±0.05 

- 

- 

- 



JJA 

0.01 ±0.06 

—2 ±2 

0±1 

- 



SON 

-0.03 ±0.03 

<1 

0±1 

- 

Ofu-olosega 

MODIS 

DJF 

-0.03 ±0.02 

1 ±2 

0±1 

- 



MAM 

-0.02 ±0.02 

3 ±2 

-1±1 

- 



JJA 

-0.01 ±0.02 

0 ± 2 

1±1 

- 



SON 

-0.01 ±0.02 

1 ±2 

0±1 

- 


AATSR 

DJF 

0.02 ±0.03 

—3 ±4 

— 1 ±4 

- 



MAM 

0.01 ±0.02 

2 ± 6 

1 ±4 

- 



JJA 

0.00 ±0.02 

0 ± 8 

0 ±4 

- 



SON 

0.01 ±0.02 

-2 ±6 

— 1 ±4 

- 

Fiji 

MODIS 

DJF 

0.04 ±0.06 

— 1 ± 3 

0 ± 3 

- 



MAM 

- 

- 

- 

- 



JJA 

0.04 ±0.06 

0 ±4 

0±1 

- 



SON 

0.04 ±0.06 

1 ± 5 

0±1 

- 


AATSR 

DJF 

-0.03 ±0.04 

— 6 ± 6 

1 ±4 

- 



MAM 

- 

- 

- 

- 



JJA 

0.06 ±0.05 

1 ± 6 

3 ± 6 

- 



SON 

0.04 ±0.05 

3 ± 6 

4 ± 6 

- 


Cloud top pressures are between 860 and 440 mb. 

1 Difference between average upwind value (—400 to 0 km) and highest (AOD) or lowest (CER) downwind value. 

2 Difference between average upwind and downwind values. 

3 “Peak distance” refers to the approximate distance of peak AOD downwind from the volcanic vent. 
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Figure 5. Scatter plots showing the difference between upwind and downwind Top of Atmosphere Short-Wave flux (a, c, e, g, i, k) for 
volcanoes and “controls”. Positive values indicate higher upward SW flux downwind of the volcano. Profiles (b, d, f, h, j, 1) indicate seasonal 
average SW fluxes for arcs of j at equal distance upwind and downwind from the volcano or island. December-February is red, March-May 
is blue, June-August is green and September-November is black. 


this decrease in CER downwind of Tristan da Cunha corre- 
lates well with a downwind increase in COD. At the active 
volcanoes (e.g. Fig. 3c or k), downwind CER remains low 
at greater distances from the vent than COD remains ele- 
vated. CER, though influenced by orographic processes, is 
more strongly affected by volcanogenic aerosol at Kilauea 


and Piton de la Foumaise, and regional trends in cloud prop- 
erties at Yasur. 
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Figure 6. Comparison of SO 2 flux from East Rift as reported by 
Elias and Sutton (2012), with AOD (550 nm) measurements from 
MODIS Aqua on the same day. The lack of relationship between 
ground- and satellite -based measurements may be due to several 
factors including dilution effects and differences in measurement 
times. Black triangles show MODIS AOD retrievals and the black 
dashed line indicates linear regression line. The red crosses show 
AOD values after de-trending for a linear relationship with wind 
speed, as demonstrated in Fig. 11. The red dashed line shows cor- 
rected linear regression. 
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Figure 7. Plots of the difference between numbers of cloud re- 
trievals (cloud fraction >0.2) downwind and upwind of volcanic 
vents (a) and island summits (b) from MODIS (Aqua), 2002-2008. 
At Hawai’i,and to some extent Reunion, there is a notable increase 
in the number of cloud retrievals made downwind relative to upwind 
within 1 00 km of the island. 


3.5 Consistency of measurements and sources of 
uncertainty 

The general trends in aerosol and cloud properties retrieved 
from MODIS and AATSR are consistent: both instruments 


show elevated AOD and suppressed CER downwind of the 
volcanoes but not control islands, and elevated COD over all 
land (Fig. 8, Table 2). However, seasonal differences vary 
between the instruments and AATSR results generally have 
much higher standard errors than MODIS due to longer re- 
peat time (3 days relative to 1 day), shorter period of cover- 
age (6 relative to 10 years) and consequently smaller sample 
size. The absolute values for CER are greater for MODIS 
than AATSR retrievals, likely due to the different penetra- 
tion depths of the instrument wavelengths. The wavelengths 
most sensitive to CER are 1 .6 and 2. 1 pm for AATSR ORAC 
and MODIS (Joint Atmosphere product) retrieval algorithms 
respectively. The shorter wavelength used by AATSR pen- 
etrates deeper into the clouds, so will sample smaller cloud 
droplets in a non-precipitating cloud, which typically have 
increasing cloud droplet size with height (Platnick, 2000; 
Sayer et al., 2011). 

Standard error in seasonal and multiannual averages is 
greatest over land for both cloud and aerosol properties (large 
error bars on Figs. 3 and 4). This is the result of (a) fewer data 
points contributing to the mean (i.e. an arc of j is shorter 
closer to the volcano) and (b) a greater spread in the retrieval 
values made over land. Observations made beyond the ex- 
tent of the island are therefore most robust. The higher un- 
certainty in retrievals made over land may also contribute to 
the “anticipation” of the volcano, seen, for example, ~ 50 km 
upwind of Kilauea in Fig. 3a and Yasur in Fig. 3b. Un- 
certainties in ECMWF wind direction will smear upwind- 
downwind differences radially in our analysis. Isolated to- 
pographic peaks, such as these volcanoes, introduce pertur- 
bations to the regional wind fields, resulting in local turbu- 
lence near the volcano itself and local differences from the 
ECMWF wind fields (spatial resolution = 1.5°) used in our 
analysis. 

As aerosol and cloud properties are both dependent on the 
same meteorological processes, there is a potential for spuri- 
ous correlation between AOD and cloud properties, unrelated 
to volcanic aerosol. For example, cyclones are associated 
with both increased cloud fraction (Field and Wood, 2007) 
and increased AOD due to high relative humidity and wind 
speeds, resulting in elevated hygroscopic growth of aerosols 
(e.g. Seinfeld and Pandis, 1998; Twohy et al., 2009) as well 
as increased sea spray emission (e.g. Grandey et al., 2011). 
Artefacts in satellite retrievals of AOD under broken cloudy 
conditions may also lead to spurious correlations between 
aerosol and cloud properties (Grandey et al., 2013), but as 
we use retrievals from the full range of atmospheric condi- 
tions such effects should not have a significant impact on our 
results. False correlations between aerosol and cloud proper- 
ties may also be introduced by analysis of satellite data over 
large regions containing significant variation in aerosol type, 
cloud regime and average synoptic conditions (Grandey and 
Stier, 2010). We limit the chance of our results being affected 
by spurious correlations due to spatial variations in climate 
by limiting measurements to regions of up to only 4° across, 
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Figure 8. Comparison of multiannual averages for AATSR (grey) and MODIS (red) data, presented as in Fig. 3 and 4. 


the largest regional size recommended by Grandey and Stier 
(2010) for analysis of AIE. Only in data from the 4° sur- 
rounding Yasur do we observe a systematic regional trend in 
aerosol properties. There were many more cloud retrievals 
(cloud fraction > 0.2, Fig. 7) downwind than upwind, though 
COD remained similar. This means that although aerosols are 
elevated downwind of Yasur, there is not evidence of statisti- 
cally significant AIE. 

4 Discussion 

4.1 Observations of volcanic aerosol 

The identification of VAIE is easiest where other sources 
of aerosol are low. We assume that for the islands dis- 
cussed here, primary sea spray aerosol dominates aerosol 
populations (e.g. Sayer et al., 2012; Smirnov et al., 2012; 
Huang et al., 2010). We test this assumption by examin- 
ing the linear relationships between mean daily AOD and 
surface wind speeds within 50 km radius of the volcanoes 
(Fig. 11). The slopes of the linear regression between AOD 
and binned wind speed take values ranging from 0.003-0.005 
for Kilauea, Yasur and La Reunion, respectively in Fig. 11. 
This is within the range found by Smirnov et al. (2012) from 


examination of ship-borne sun photometer AOD measure- 
ments and near-surface wind speeds (0.004-0.005). This sup- 
ports the assumption that sea spray is the most important 
background aerosol at these volcanoes, rather than, for ex- 
ample, organic carbon or mineral dust from the islands. 

Time series of daily AOD from MODIS data all show a 
large day-to-day variability. This is partially attributable to 
the amount of both sea spray and volcanogenic aerosol be- 
ing heavily dependent on synoptic conditions. Although the 
emission of volcanic gases is independent of atmospheric 
conditions, the oxidising agents and reaction rates for the 
formation of SO 4 - from SO 2 are dependent on atmospheric 
conditions such as supersaturation and solar irradiation lev- 
els (Eatough et al., 1995). The formation of aerosol may take 
place in a different manner in clear and cloudy conditions, 
and dilution effects will also vary with wind speed and di- 
rection. The emission of volcanic gases also varies through 
time. While at Kilauea, SO 2 emission is passive and accom- 
panies long-lived flows and lava lake activity, at Yasur gas 
flux and composition varies more rapidly according the stage 
of Strombolian explosion. At Piton de la Fournaise inter- 
eruptive flux is low, but SO 2 emission accompanies frequent 
minor explosions and lava flows. There is a general corre- 
lation between net downwind aerosol at the three volcanoes 
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Figure 9. Scatterplots show the average of rotated MODIS Aqua 
AOD (a, b) and CER (c, d) at Piton de la Fournaise during inter- 
eruptive (a, c) and eruptive (b, d) periods, 2002-2008. Cloud top 
pressures are 860-440 mb. The volcano summit is at (0,0) and ar- 
rows indicate wind direction. Downwind is the upper quadrant in 
each case, n values refer to the average sample size in each image. 
Only data where sample sizes exceed 1 5 and relative standard devi- 
ations are less than 10 % are shown. 


and average background (inter-eruptive) SO 2 flux (e.g. Ta- 
bles 1 and 2). 

4.2 Volcanic aerosol and cloud interaction 

The impact of volcanogenic aerosol depends on synoptic 
conditions and the state of any pre-existing cloud. Secondary 
sulfate aerosol formed from ambient reactions of volcanic 
SO 2 are typically <0.1 pm in diameter, and are therefore too 
small to act as CCN where the level of supersaturation is low 
(e.g. Pierce and Adams, 2007). Aerosol from shipping most 
commonly results in tracks of lower CER where the aerosol 
plume is able to mix directly with an overlying cloud deck 
before major dispersal. Typical ship tracks form when the 
boundary layer is shallow (< 800 m), wind speeds are mod- 
erate, relative humidity is high and the difference in temper- 
ature between air and sea is low (Durkee et al., 2000). The 
impact of additional CCN on orographic cloud has been less 
well studied, and conditions most likely to result in cloud 
brightening are unclear (Muhlbauer et al., 2010). 

All of the islands investigated here are in regions where 
the free atmosphere is dominated by trade winds, except Tris- 
tan da Cunha, where westerlies dominate. We expect that the 
measurements of cloud properties are most commonly from 
decks of marine stratocumulus over the oceans, with contri- 
butions from orographic cloud over land and in the islands’ 


wakes. As our results consist of seasonal and multiannual 
averages of retrievals, they contain contributions from days 
with a range of meteorological conditions (Mean COD and 
cloud top pressures for MODIS data 2002-2013 are shown 
in the Supplement). The mean values for COD and CER pre- 
sented here therefore capture net conditions and may not bear 
a resemblance to the atmospheric processes on any particular 
day. Our measurements of time-averaged impact would be 
complemented by case-by-case analysis of volcanic tracks 
(e.g. Gasso, 2008) to improve our understanding of the influ- 
ence of atmospheric conditions and emission characteristics 
on cloud properties. 

4.3 Evidence for volcanic aerosol indirect effects 

The anticorrelation between AOD and CER values down- 
wind of the degassing volcanoes provides strong evidence 
for a first indirect effect at Kilauea and Piton de la Fournaise. 
This is confirmed by the perturbation to ToA upward SW 
flux downwind of the active volcanoes seen in CERES data. 
At all the active volcanoes minimum CER is found at similar 
downwind distance as peak AOD in averaged retrieval data. 
The difference between background (upwind) CER and the 
minimum value found downwind is greatest at Kilauea (up 
to —8 pm, decrease of 35 %) and lower at Yasur and Piton 
de la Fournaise (up to —4 pm, decrease of 18%). This dif- 
ference far exceeds the expected uncertainty in MODIS CER 
retrievals of 2-3.5 % (Platnick et al., 2003). In common with 
the persistent downwind elevation in AOD, CER at Kilauea 
and Yasur remains lower than upwind values at least as far as 
400 km downwind, while CER returns to its average upwind 
value within 100 km at Reunion. Droplets in “clean” clouds 
upwind of the volcanoes have average droplet effective radii 
of 21-22 pm. 

The presence of secondary AIE are harder to demonstrate. 
Liquid water path (LWP, kgm -2 ) reaches its peak value 
within ~ 50 km of the volcanic vents and is very slightly el- 
evated downwind relative to upwind at the volcanoes (Fig. 3 
and 4). This could indicate an increase in cloud lifetime due 
to drizzle suppression (e.g. Lohmann et al., 2005), but it 
may also be the consequence of the contribution of cooling 
and condensing water vapour emitted from the volcano. At 
Kilauea, and to some extent Piton de la Fournaise, there were 
an elevated number of cloud retrievals with cloud fraction 
>0.2 downwind relative to upwind (Fig. 7). Water vapour 
flux at Yasur exceeds 13 x 10 3 Mgday _1 (Metrich et al., 
2011). Emissions of water vapour into the atmosphere at 
Kilauea and Piton de la Fournaise will be even greater due 
to the additional contribution of evaporated seawater where 
lava flows meet the sea (e.g. Edmonds and Gerlach, 2006; 
Gouhier and Coppola, 2011). The condensation of evapo- 
rated seawater may result in an increase in LWP content 
downwind that would mask any evidence of secondary AIE. 

In spite of the strong correlation, our MODIS and AATSR 
measurements do not allow a direct measurement of the 
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Figure 10. Profiles showing average COD (a-c) and CER (d-f) at Kilauea, Yasur and Piton de la Fournaise, separated according to Froude 
number squared {Fr 2 ) as estimated from ECMWF data. Fr 2 = 0-0.1 is blue, Fr 2 = 0. 1-0.5 is green and Fr 2 >0.5 is red. 



wind speed at sea level (m/s) 


Figure 11. Daily average AOD at 550 nm measured from MODIS, 
Aqua between 2002 and 2008 with 50 km of the coasts of (a) 
Hawai’i, (b) Yasur and (c) Reunion plotted as a function of hori- 
zontal wind speed at sea level from ECMWF. Error bars show the 
standard deviation in AODs for bin intervals of 1 m s _1 . 


Twomey effect because retrievals of aerosol and cloud prop- 
erties are mutually exclusive (i.e. aerosol properties are re- 
trieved only where pixels are not flagged as cloud). However, 


the increase in average downwind ToA upward SW flux seen 
in CERES data for all sky conditions shows a radiative im- 
pact over a similar area to the downwind cloud and aerosol 
perturbations, providing additional evidence of a first indirect 
effect associated with volcanic aerosol. 

5 Extrapolating to global volcanic aerosol indirect 
effects 

Elevated AOD and suppressed CER are observed over a 
decade downwind of three persistently active volcanoes 
with different eruptive characteristics: Kilauea (strong de- 
gassing), Yasur (Strombolian eruptions) and Piton de la Four- 
naise (minor explosions and lava flows). Top of Atmosphere 
Short Wave radiative flux is apparently elevated by at least 
10 Wm' 2 downwind of all three degassing volcanoes at dis- 
tances >150 km, and is even higher within 100 km downwind 
of Kilauea and Piton de la Fournaise. Time-averaged cloud 
data at three “control” islands showed no significant suppres- 
sion of CER downwind of the islands or perturbation to ToA 
upward SW flux exceeding ±10 Wm -2 . 

Our approach builds on previous studies of volcano tracks 
during periods of elevated activity (Gasso, 2008; Yuan et al., 
2011) by estimating average volcanic impact over 6-10 
years. This is a step towards the measurement of the long- 
term impact of persistent volcanic activity in the present- 
day atmosphere. In the three test cases described here, we 
observe some general trends that could be extrapolated to 
other volcanoes. High average downwind AOD is linked 
with high background SO 2 flux. Although daily ground- 
based measurements of SO 2 do not necessarily correlate with 
daily satellite retrievals of AOD, average background SO 2 
flux at the three volcanoes is proportional to time-averaged 
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upwind-downwind difference in AOD (Tables 1 and 2). Both 
the greatest upwind-downwind differences in AOD and CER 
and the greatest distance to minimum CER were found at 
Kilauea, the volcano with the greatest SO 2 flux. 

Of the 49 continuously erupting volcanoes for which An- 
dres and Kasgnoc (1998) present SO 2 fluxes, only 8 are is- 
lands further than ^50 km from the nearest land mass and to- 
gether contribute ~ 14 % of the SO 2 flux in the compilation. 
Our observations at Kilauea and Yasur alone capture the local 
impact of about half of these emissions. The remaining SO 2 
flux from continuously degassing volcanoes originates on the 
continents (24 volcanoes) or within ^50 km of the nearest 
large landmass (17 volcanoes). These “coastal” volcanoes, 
including Etna, Bagana and Sakura-jima, emit almost ~ 50 % 
of global continuous emissions (Andres and Kasgnoc, 1998). 
Observations of AIE in continental or coastal settings may be 
achievable using trajectory analysis of air parcels from the 
volcano to allow comparison of aerosol-laden atmosphere to 
average background conditions. Our observations of VAIE at 
Piton de la Foumaise demonstrate that active volcanoes with- 
out a constant SO 2 flux to the atmosphere also have an im- 
pact on time-averaged local cloud properties. A complete in- 
ventory of tropospheric cloud alteration by volcanoes should 
therefore also include emissions associated with minor ex- 
plosions. 

The VAIE observed in this study is unlikely to be represen- 
tative of the average present-day impact of volcanoes world- 
wide. Importantly, the impact of additional aerosols on cloud 
microphysics is greater for pristine clouds than for polluted 
regions (e.g. Lohmann et al., 2005; Rosenfeld et al., 2008). 
The change in cloud droplet size downwind of Kilauea, Ya- 
sur and Piton de la Foumaise is expected to be large rela- 
tive to that for volcanoes that emit aerosol into heavily pol- 
luted cloud decks (e.g. Etna, Masaya). However, these obser- 
vations of the most pristine parts of the Earth’s present-day 
atmosphere are likely to be the best observable analogue of 
pre-industrial aerosol effects. 

The Supplement related to this article is available online 
at doi:10.5194/acp-14-10601-2014-supplement. 
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